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■ Abstract 

The goal of this article is to provide a practical method to calculate, in a scalar 
pi i' theory, accurate numerical values of the renormalized quantities which could 

be used to test any kind of approximate calculation. We use finite trunca- 
tions of the Fourier transform of the recursion formula for Dyson's hierarchical 
model in the symmetric phase to perform high-precision calculations of the 
unsubtracted Green's functions at zero momentum in dimension 3, 4, and 5. 
I We use the well-known correspondence between statistical mechanics and field 

theory in which the large cut-off limit is obtained by letting (3 reach a critical 
value (3c (with up to 16 significant digits in our actual calculations). We show 
that the round-off errors on the magnetic susceptibility grow like {j3c — (3)~^ 
near criticality. We show that the systematic errors (finite truncations and 
volume) can be controlled with an exponential precision and reduced to a 
level lower than the numerical errors. We justify the use of the truncation 
for calculations of the high-temperature expansion. We calculate the dimen- 
sionless renormalized coupling constant corresponding to the 4-point function 
and show that when /? — > /?c, this quantity tends to a fixed value which can 
be determined accurately when D = 3 (hyperscaling holds), and goes to zero 
like {Ln{f3c - f3))'^ when D = 4. 
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I. INTRODUCTION 



Finding closed form, exact analytical solutions to difficult problems is considered as a 
great achievement in theoretical physics. In recent years, the development of fast computers 
and of easy electronic communications has enlarged the class of solutions which can be 
considered as completely satisfactory. Lengthy expressions can be manipulated symbolically 
or numerically and communicated to others (for a concrete example, see for instance the 
Appendices of Ref. [0). Sometimes, the solution of a problem requires a combination of 
iterations and expansions which can be performed with any desirable precision in a short 
amount of time using a friendly environment such as Mathematica. An example of such 
a solution is the calculation of the spectrum of the one-dimensional quantum anharmonic 
oscillator described in Ref. 0. In this example, even though no closed expression for the 
eigenvalues and eigenfunctions is available, the beginning of the spectrum can be obtained 
numerically with great precision and almost instantly using Mathematica. We could say that 
this problem is numerically solvable. 

Scalar field theory on a Euclidean lattice is a difficult problem with many important 
applications, such as the interactions of strongly interacting light particles (pions, kaons,...), 
the generation of mass in the standard model of elementary particles, and the theory of 
critical phenomena. The renormalization group method helps us to understand the nature 
of the continuum limit 0] for such a model. However, in the case of short range (nearest 
neighbor) interactions, we are still far away from the numerical solvability mentioned above. 

In order to solve scalar field theory, one needs an approximation scheme such that: a) 
the zeroth-order approximation preserves the main qualitative features of the model, b) 
the zeroth-order approximation is analytically or numerically solvable, and c) the zeroth- 
order approximation can be improved systematically and in a practically implementable 
way. The fact that Wilson's approximate recursion formula satisfies the requirement a) is 
justified in Ref. Q. The approximate recursion formula is an integral equation with one 
variable which can be handled by standard numerical methods. The main sources of errors 
are the finite number of points of integration and the parametrization of the behavior of 
the tails of the functions integrated. The errors can be reduced by reducing the size of 
the tails and increasing the number of points. We found this approach time consuming 
and inefficient. However, using the Fourier transform of the recursion formula, we found 
a natural approximate method with a fast implementation and a control of the systematic 
errors which is better than exponential. This method is the main calculational tool used 
and discussed in the present article. 

The approximate recursion formula is closely related to the recursion formula appearing 
in Dyson's hierarchical model More precisely at fixed dimension, there exists a 
one parameter family of recursion formulas which interpolates continuously between the 
two. Following Ref. |0, we call this parameter (. Seen from a practical point of view 
in our calculation, the extension to another value of ( in the recursion formula amounts 
to changing one number in one line of a two page program. Physical quantities such as 
the critical exponents vary slowly when ( is varied, but there is nothing that singles out a 
particular value of (. In the following, we specialize the discussion and the numerical study 
to the case of Dyson's model = 1/D) because this model has been studied [[7|-[l0[1 in great 
detail in the past. Dyson's model is a well-defined lattice model and it admits the same 
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kind of expansions (weak and strong coupling, large- etc.) as any other scalar model on a 
cubic lattice. The basic recursion formula and the approximation methods are explained in 
section 0. 

For simplicity, all the calculations done below use an initial Ising measure. The infinite 
cut-off limit is obtained by fine-tuning the only adjustable parameter, namely, the inverse 
temperature (3. All the results will be expressed in terms of (3c — (3. The general procedure 
which relates small values oi (3c — (3 to large values of the UV cut-off is well-known 
and will not be repeated here. From a calculational point of view, the discussion would be 
essentially the same if instead we had fixed (3 = 1 and considered initial measures depending 
on a cut-off and several bare parameters. 

The rest of the presentation is based on the following empirical fact: the systematic 
errors due to finite volume fall exponentially with the number of iterations {nmax and the 
systematic errors due to finite dimensional truncation fall faster than exponentially with 
the dimension of the truncated space {Imax)- In calculations involving double precision, one 
can, without encountering major difficulties, choose rimax and l^ax in such a way that if 
we increase these parameters further, no change is observed in the results. One can then 
first determine the critical temperature and the numerical errors. We need to discuss these 
first, because the numerical errors are a fundamental limitation whenever the precision of 
the arithmetic operations is fixed, and there is no point in trying to reduce the systematic 
errors much below the numerical errors. 

The phase structure of the approximated models is discussed in section |T|. We show 
that, for understandable reasons, there is nothing that can be identified as a low-temperature 
phase. We also give a practical method to identify the inverse of the critical temperature, 
denoted /3c, with an optimal accuracy. In the following sections, it will always be assumed 
that we work in the symmetric phase {(3 < (3c)- 



The numerical errors on the magnetic susceptibility are studied in section where we 
show that the relative errors obey the approximate law 

where 6 is the precision used to perform the arithmetical operations. We also show that this 
law follows from an approximate renormalization group calculation. As a by-product we 
obtain a numerical estimate of the critical exponent 7 in good agreement with the existing 
estimates for D = 3. As a general remark, we have a much better quantitative control 



on the details of renormalization group arguments in D = 3 than in D = A, where confluent 
logarithmic singularities make the analysis more delicate. 

We then discuss the systematic errors. The volume dependence of the magnetic suscep- 
tibility is discussed in section 0. We show that in order to calculate the susceptibility at an 
inverse temperature (3, in D dimensions, and with relative errors A, one needs a number of 
lattice sites which is of the order of (A(/9c — (3)^)~^^'^ . We then show that when the Fourier 
transform of the recursion formula of Dyson's hierarchical model is projected into a finite 
dimensional space of dimension Imaxi the relative errors on the susceptibility associated with 
this truncation decrease faster than e~" for some positive number a. We show in sec- 
tion |V11| that similar results apply to the high-temperature expansion of the susceptibility, 



justifying the procedure used in Ref. [T^ • 
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At this point, we have only discussed the susceptibihty, or in other words, the 2-point 
function. The 4-point function is also provided by the calculational method at no extra 
cost. However, the calculation of the corresponding renormalized coupling constant requires 
a subtraction. In the discussion of the errors on the susceptibility, we have explained that 
while we are iterating the recursion formula we lose significant digits "from the right" . With 
the subtraction, we also lose significant digits "from the left". This is explained in section 
Vlll|, where we calculate a "dimensionless coupling constant" inspired by the field theory 



definition of Ref. ||T2[ and denoted A4 . 

In the case D = 3, we show that when (3 — > /5c, A4 tends to a fixed non-zero value that 
we were able to calculate with 6 significant digits. This is very convincing evidence that 
hyperscaling holds for the model considered here. For the nearest-neighbor Ising model on 
a 3-dimensional cubic lattice, it is very hard to decide if hyperscaling holds on the basis of 
Monte-Carlo simulations [0 or high-temperature expansion [|14|. In the present case, it is 
a short and straightforward calculation. This shows that it would be worth trying to in- 
terpolate perturbatively between Dyson's model and nearest neighbor models, transforming 
the qualitative approximation into a quantitative approximation. The accuracy of A4 
decreases when D increases. However, we were able to obtain good evidence that in D = 4, 
A4 decrease like (Ln(/3c — in good agreement with the behavior obtained with the field 
theory method at lowest order in perturbation theory. 

Our results show that it is possible to calculate accurately, and without major effort, the 
renormalized quantities which can be extracted from the 2- and 4-point functions. These 
calculations can be done for an extended range of small values of /3c — (3- In other words, for 
existing computers, requirement b) is fulfilled by the hierarchical approximation provided 
that one does not require too-small values of (3^ — (3- Numerical solvability can sometimes 
completely change our point of view regarding a problem. Taking the example of differential 
equations, at a time when their numerical solutions were not achievable, one could dream 
of a perfectly deterministic approach to natural phenomena in which everything could be 
known once the rules of evolution and the initial conditions were given. 

This paper provides a calculational method to anyone who would like to check an approx- 
imation method with accurate numerical results. As explained above, most of the approx- 
imation schemes which apply to a scalar field theory with nearest neighbor interactions on 
a cubic lattice also apply to Dyson's model. We are presently testing the validity of several 
well-known but not easy-to-control expansions: the renormalized perturbative expansion, 
the loop expansion, and the large- expansion. 



II. THE RECURSION FORMULA AND ITS FINITE DIMENSIONAL 

TRUNCATIONS 

In this section, we briefly describe Dyson's hierarchical model and the methods used to 
calculate the average of arbitrary powers of the zero momentum component of the scalar 
field. The models considered here have 2^"^'^'^ sites. We label the sites with Umax indices 

Xn^ax each index being or 1. In order to visualize the meaning of this notation, 

one can divide the 2"™""' sites into two blocks, each containing 2"™""="^ sites. If Xn^^x = 
the site is in the first block, if = 1, the site is in the second block. Repeating this 

procedure Umax times (for the two blocks, their respective two sub-blocks, etc.), we obtain 
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an unambiguous labeling for each of the sites. Two sites differing only by xi are in the same 
block of size 2. We write the action as 



1 T^max 

^ = -2 ^(4)" E ( E ■ (2.1) 

n — 1 ^rimax vi^n + l i^l 

The index n, referred to as the "level of interaction" hereafter, corresponds to interactions 
of the total field in blocks of size 2'^. The free parameter c, controlling the strength of the 
interactions, is set equal to 2^"^/^ in order to approximate a nearest neighbor model in 
D-dimensions. In this article, we will consider the cases D = 3, 4, and 5. 

The field is integrated with a local measure PFo(0) which needs to be spec- 

ified. In the following, we use an Ising measure, VFo(0) = ^{(f^ — 1), which takes only the 
values ±1. To the best of our knowledge, the results presented do not depend on this specific 
choice and would also apply, for instance to the case of Landau-Ginzburg measures of the 
form Wo{(t)) = e-^'^^-^'^\ 

The integrations can be performed iteratively using the recursion formula 



where C„+i is a normalization factor which can be fixed at our convenience. The relation 
between this recursion formula and Wilson's approximate recursion formula is discussed 
in Ref. 0. Introducing the Fourier representation 

Wn{<P) = I %'"^W^{k) , (2.3) 
and a rescaling of the source by a factor 1/s at each iteration, through the redefinition 

Rnik) = Wr,i\) , (2.4) 

s 

the recursion formula becomes ^ 

Rn+iik) = ^M-l/3ils'r^'-^)iRni-jr . (2.5) 

The rescaling operation commutes with iterative integrations and the rescaling factor s can 
be fixed at our convenience. 

In the following, we fix the normalization constant C„ in such way that -Rn(O) = 1. Rn{k) 
has then a direct probabilistic interpretation. If we call M„ the total field J2 4'x inside blocks 
of side 2", and < ... >„ the average calculated without taking into account the interactions 
of level strictly larger than n (or in other words, as if n were equal to Umax), we can write 

We see that the Fourier transform of the local measure obtained after n iterations generates 
the zero-momentum Green's functions calculated with 2" sites, and can thus be used to 
calculate the renormalized coupling constants at zero momentum. 
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The choice of s is a matter of convenience. For calculations in the symmetric (high- 
temperature) phase not too close to the critical temperature, or for high-temperature ex- 
pansions PJ^, the choice s = a/2 works well. For calculations very close to the critical 
temperature, the choice s = 2c~2 prevents the appearance of very large numbers. In the 
following calculations, we always use s = 2c~2. 

In the following, the finite volume magnetic susceptibility is defined as 

UP) ^ ^^f^ . (2.7) 
With the rescaling s = 2c^^, we have the relation 

Xn = -2an,i{'^r ■ (2.8) 

The initial condition for the Ising measure is Rq = cos{k). For the Landau- Ginsburg mea- 
sure, the coefficients in the fc— expansion need to be evaluated numerically. 

In previous publications @Jl3], we used the P expansion of Eq. ( p.5|) to calculate the 



high-temperature expansion of the magnetic susceptibility of Dyson's hierarchical model up 
to order 800 with an Ising or a Landau- Ginzburg measure. The calculation of the large- 
order coefficients requires a lot of computing time. We found that using a truncation in 
the expansion in at order 50 could cut the computer time by a factor of order 100 while 
having effects on the values of the coefficients which were smaller than the errors due to 
numerical round-off. 

We consider the finite dimensional approximations of degree Imax'- 

Rn{k) = 1 + an,ik^ + an,2k^ + + a„,Ua.^''™^ • (2-9) 

After each iteration, non-zero coefficients of higher order {an+i,imax+i ^tc.) are obtained, but 
not taken into account (i.e. set to zero as part of the approximation) in the next iteration. 
More explicitly, the recursion formula for the a„ „ reads: 



^lm,ax(Y^ „ „ N (2/)! (c\l( 1 n\l-m 

^n+l,m — „; n C9/V / / TTTl ' \Z,-l\J) 



III. THE PHASE STRUCTURE OF THE APPROXIMATED MODELS 

From a field theoretic point of view, an important feature of the hierarchical model is its 
second order phase transition. Our first task will be to identify Pc for the well-studied PJT6[| 
case of an Ising measure. The truncation described above can be used to calculate ^ exactly 
the high-temperature expansion of the magnetic susceptibility up to order l^ax — 1, and one 
can think about the truncation as a partially re-summed high-temperature expansion. It 
seems thus unlikely that when Imax becomes large, we would obtain sensible results when 
(3 > (3c- This guess is confirmed by empirical numerical experiments at fixed (3. 

The truncated recursion formula shows very clearly the existence of a high-temperature 
phase where for j3 smaller than Pc and n large enough, we have the scaling law 
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fln.l OC 2 o . 



(3.1) 



Given the choice of scahng factor s = 2c 2 and the definition of c discussed in section ||, 
this scahng is equivalent to the "central limit" behavior 



< (M^)'" >„oc 2"^ . 



(3.2) 



_2_ 

' D , 



This situation is characterized by ratios a„+i,i/a„^i reaching the value c/2 = 2 

On the other hand, when f] is increased sufficiently, there is a sudden change of behavior. 
However, there is nothing that we can identify with a low-temperature phase |T2|. In all 
the cases that we have considered, the change is signaled by the fact that when n increases, 
the ratio an+i,i/0'n,i becomes larger than 1, and then starts making unpredictable changes, 
never returning to any kind of behavior like Eq. ( |3.2|) or the behavior characteristic of the 
low temperature phase, namely 



< {Mr, 



^2q 



>n.0(. 



(3.3) 



The "irreversibility" of this process allows us to identify unambiguously /?c, in the sense that 
a calculation at finite n gives upper and lower bounds on /3c. By increasing n, we can obtain 
sharper bounds. This procedure is illustrated for D = 3 in Fig. |1|. 

We see that a calculation for n up to 50 allows us to resolve the 10-th digit of /3c, and a 
calculation for n up to 60 allows us to resolve the 11-th digit. Proceeding similarly, we can 
determine the numerical value of with as many significant digits as the computational 
method allows. Double precision Fortran calculations made with Imax = 80 are reported 
below for the Ising model in 3, 4, and 5 dimensions. The results are in agreement with the 
bounds found in ref. |jl6[ with independent (and exact) calculational methods. The third 



column gives the minimal value of n which allows a resolution of the 16 significant digits. 



D 



TTKLXfYLin 



3 
4 
5 



1.179030170446270 
0.6654955715318593 
0.4569633006170210 



102 
111 
132 



32 
43 
45 



(3.4) 

Subsequently, l^ax was lowered by small steps until the value of Pc changed. This exper- 
iment shows that the change occurs at values much smaller than 80. In the fourth column of 
Table 1, we give the minimal value of l^ax such that the stable value of (3c with 16 significant 
digits can be reached. 

One may wonder if the precise value of f3c is dependent on the numerical aspects of 
the calculation such as the round-off errors, which will be discussed in the next section. 
To settle this question, we have used methods which perform the arithmetic operations 
differently (these methods are explained in detail in the next section) and found the very 
same values of (3c- In conclusion, we have found a reasonably robust value of (3c which is 
consistent with existing results. 

We have calculated /5~^ for dimensions much larger than 4. The results are displayed in 
Fig. 1^. We see that this quantity grows linearly with the dimension. The explicit calculation 
of the high-temperature expansion |jl| suggests that when D 00, i.e. when c — * 2, we have 



7 



- ' (4-c)(2-c) )'° • P-^' 

which imphes that 

This estimate of the slope is in good agreement with the data. We found a slope of 0.714, 
while {2Ln{2))~^ = 0.721. Next we will study the various sources of error occurring when 
one approaches f3c from below. 



IV. THE ROUND-OFF ERRORS 

Round-off errors can play an important role when recursive methods are used, because 
they may grow faster than the improvement of the results due to the repeated use of the 
method. For this reason, we have studied them with three independent methods. By 
independent, we mean that the arithmetic is performed in a completely uncorrelated fashion. 

We have compared our original Fortran calculation on a DEC-alpha with three other 
calculations. The first one was the same program run on a MIPS. The second one was a 
Mathematica program where a higher precision in the arithmetic operations was set. The 
precision was adjusted in such a way that the susceptibility was obtained correctly with 16 
significant digits. Thirdly, we have compared the calculation with the one obtained with a 
slightly different rescaling, namely s = 1.98c~2, a method already used in Ref. P,p^. All 
these calculations were performed with l^ax = 60, which is beyond what we need (see next 
section). 

The relative differences in the finite volume susceptibility are shown in Fig. 0, for D = 3 
and Pc — P = 10~^^. The figure shows clearly that the three types of discrepancies are 
essentially the same. Since the three types of errors are uncorrelated, we can identify them 
with the round-off errors and calculate them with the most convenient method. Using the 
third method, we have calculated the round-off errors for various values of /3. In all the 
cases, the logarithm of the relative error grows linearly at the beginning and then stabilizes 
at a constant value. The period of linear growth corresponds roughly to the iterations where 

2 

Xn+i — 2"DXn- For larger n, the value of the susceptibility stabilizes, with changes decreasing 

2 

by a factor 2~d at each step. During this second stage, the numerical errors do not grow 
significantly. 

We now proceed to discuss the asymptotic values of the errors, in other words, the 
stable value they reach for n sufficiently large. We have collected these values for various 
temperatures and D =3, 4, and 5 in Fig. This shows that the relative error is in good 
approximation 10~^^(/5c — independently of D. 

These empirical results have a simple explanation in terms of the linearized theory. 
Suppose that 5 is a typical round-off error in a calculation (e.g. 10~^^), and that A is the 
largest eigenvalue of the linearized renormalization group transformation near a given fixed 
point. One expects the numerical error on a„ 1 to be of the order A'^^. With the rescaling 
used in this paper, this means that the errors on the susceptibility are of the order 

I^Xnl ~ A"5(-)" . (4.1) 



8 



(4.2) 



Now for n such that A" ~ {j3c — P) ^, we have x ~ (~)"- Plugging this into the previous 
equation, we get 

I — I ^ 

which is the empirical result found above. 

For D = 3, one can check the details of the above argument and, as a by-product, obtain 
an estimate of A. The slope of the increasing part of Fig. |^ is approximately 0.154, with an 
estimated error of order 0.001. This implies a value 10°-^^^ = 1.427 for A. Using the usual 
formula for the critical exponent, 

(A 



we obtain the value 7 = 1.30, in good agreement with existing estimates [p|-[TT[|. 

For D = 4, the same procedure gives an exponent which is too high by about 3 per- 
cent (compared to the trivial value). This is a typical error when one does not take into 
account the marginal direction and the (related) confluent logarithmic singularity in the 
susceptibility. 



V. VOLUME EFFECTS 

In the two previous sections, we developed a qualitative understanding regarding the 
finite volume susceptibility Xn, or in other words regarding the way the susceptibility depends 
on the number of iterations. Volume effects can be important in the determination of the 



critical exponents. For instance, in Ref. |]I6[, exact calculation with almost a million sites 
gave errors of more than 10 percent in the exponent 7. We are now ready to get a better 
quantitative understanding of these effects. 

If we consider the evolution of XniP) when n increases, with (3 fixed slightly below (3c, 
we see from Fig. |1| that when we are close to the fixed point, Xn+i — '^^Xn- This lasts 
until the right order of magnitude ( ~ (/3c — ) is reached. For larger n, the value of 
the susceptibility stabilizes, with errors decreasing at each step. In this second regime, the 
measure becomes asymptotically Gaussian, and one can estimate the change in Xn from the 
change in the k"^ term. From the basic formula (2.5), one gets the estimate for the relative 
change: 

A I Xn+l Xn I o— — n /r i \ 

n = I 1 OC 2 D . (5.1) 

Xn 

From these considerations, we find the number n{i3,A) of iterations necessary to calculate 
the susceptibility at fixed /3, with a relative precision A (defined as in Eq. ( [5.1D ) : 

n{P,A) = -{^^£^){Log,oiA)+^Log,o{Pc-P)) • (5.2) 

The comparison with a numerical calculation where we required A = 10~^^ is given in Fig. 
^ for D =3, 4, and 5. The agreement with the estimate of Eq. ( |5.2|) with 7 = 1.3 (1.0) for 
D = 3 {4 and 5), is quite good. 
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The fact that we were able to stabihze sixteen digits of the susceptibihty does not mean 
that the results have sixteen digit accuracy. The asymptotic stability of the numerical 
results comes from the fact that the r. h. s. of Eq. ( |5.1j ) will go to zero whenever Xn quits 
growing. This occurs independently of the fact that numerical errors may occur while Xn 
builds up its bulk value. Consequently, a more realistic approach would be to require a 
precision consistent with the round-off errors discussed in the previous section. Imposing a 
temperature-dependent requirement A = 10~^^{j3c — we obtain values of n shown in 

Fig. |. 



VI. CONTROLLING THE EFFECTS OF FINITE DIMENSIONAL 

TRUNCATIONS 

In this section, we study the Imax dependence of the magnetic susceptibility for P < 
Pc- For notational purposes, we call x^''* the susceptibility corresponding to a given value 
Imax = I- For each calculation, the value of Umax has been increased until no change could 
be observed. The results are displayed in Fig. ^ for /?c — = 10~^. For low /, x^'^ grows 
at a not- very regular rate and within the bounds 1 < x^^^^Vx^^^ < 10. When / gets close to 
20, x'-'-' starts stabilizing with a precision which seems to be exponential. For instance, for 
= 3, the relative errors fall approximately like lO""'^'. This exponential rate is based on 
the assumption that the logarithm of the relative errors falls linearly. However, a closer look 
shows that it falls slightly faster. This is illustrated in Fig. ^ The best parametrization 
that we have found is a linear function times the logarithm of /. A more detailed analysis 
shows that this new parametrization reduces the square root of the sum of the square of 
the relative differences ((fit-data) /data) by one order of magnitude. This suggests that one 
should try to derive rigorous bounds where the errors are proportional to some inverse power 

of (^max'')' 

We have thus studied the logarithm of the relative differences (due to the change in Imax) 
divided by the logarithm of Imax for various temperatures with D = 3. The results are shown 
in Fig. 7. We then used linear fits for the part falling linearly. In other words, we assumed 
the approximate law 

. (6.1) 

X 

The results can be summarized as follows. The slope is almost independent of (3 and takes 
the approximate value -0.41 with changes of order 0.01. The intercept grows linearly with 
—Logio{Pc — P), as shown in Fig. |^. A linear fit of this data gives an intercept of the form 
1.7 — 0.83 X LogioiPc — P)- If we neglect the slow logarithmic variations and approximate it 
by a constant central value in the falling part of Fig. |^, we obtain the approximate law 

1 — 1 ~ 3.2 X 10^ X (/?e - P)~^-^ X (4.1)-'— . (6.2) 
X 

If we require these errors to be smaller than the numerical errors, we find that l^ax = 40 
is a safe choice for all the values of Pc — P accessible with double precision. Slightly larger 
values are obtained for D = 4 and 5, which confirms that the last column of Eq. ( p.4| ) 
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represents approximately the values of Imax above which no significant changes are observed. 
In conclusion, for calculations using double precision, the choice Imax = 50 is convenient and 
safe for the three values of D considered above. 

Having an acceptable control on the susceptibility guarantees that we have an acceptable 
control on the higher moments, < M^'' >„ /2'^^ for g > 1, since to leading order in the 
volume, these quantities are dominated by the disconnected parts. The precision which 
can be achieved on the connected parts (which enter in the definition of the renormalized 



coupling constants) is a more delicate question, which is discussed in section |VIII 



VII. EFFECTS OF FINITE DIMENSIONAL TRUNCATIONS ON THE HT 

COEFFICIENTS 



In a previous publication [|^, we used the truncated algorithm to calculate 800 coeffi- 
cients of the high-temperature expansion of the magnetic susceptibility. We claimed that 
this truncation did not affect the numerical values obtained. In this section, we provide a 
more systematic justification of this procedure. 

We examine the /^ax-dependence of the high-temperature coefficients of the suscepti- 
bility, for dimensions 3, 4, and 5. As in section |V^, we replace Imax with / for notational 
purposes. We denote the high-temperature coefficients as bm for the mth coefficient. 

For / and m large enough, we find good linear fits in I for the quantity 

Lniibm - y^/hm) 



Ln{l) 

where is the truncated version of the exact bm- Fig. ^ show these lines in D = 3, for 
m = 200, 300, and 400. Very similar numerical values are obtained for D = 4 and D = 5. 
For this reason, it was impossible to display the values for the three chosen dimensions in a 
single graph. The graphs in D = 4 and D = 5 are similar looking and show a linear behavior 
as good as in D = 3. We can thus express the truncated coefficients as 

6« = 6^(l-rW'+''), (7.1) 

where s and q are, respectively, the slope and intercept of the corresponding fitted line. For 
the three chosen dimensions, the lines seem to "focus" in one point close to the / = axis. 
The intercepts are approximately independent of m and take the approximate values 3.4 
{D = 3), 2.3 {D = 4), and 1.7 (-D = 5). For the slopes, we find straight line fits if we plot 

Ln{—s{m)) 
Ln{m) 

versus Ln{m) for D = 3, and versus m for D = 4 and 5. In Figs. |TT] and |T2] we have used 
every tenth coefficient in the range from m = 300. ..400. From these fits, we find for D = 3: 

S = _^.013LnM-.22_ (73) 

For D = 4 we found 
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s = _rn-""'-'\ (7.3) 

and for D = 5: 

In D = 3, for example, these results show us that we need / larger than 34 for the error 
on feiooo to be less than 10~^^, while for D = 4, / > 38. Therefore, the value of / = 50 we 
have previously used in calculating the first 800 coefficients is more than adequate. 



VIII. CALCULATION OF SUBTRACTED QUANTITIES 

In this section, we discuss the numerical aspects of the calculation of subtracted quan- 
tities. We specialize the discussion to the calculation of the "dimensionless renormalized 
coupling constant" |jl2| corresponding to the 4-point function. 



From Eq. ( |2.6| ), it is clear that the calculation of Rn{k) allows us to determine the renor- 
malized coupling constants. The first step in the calculation of these quantities is to extract 
the connected parts. In other words, we first subtract the disconnected parts from the 2k- 
point function. From a numerical point of view, this is not a trivial operation, because the 
subtracted quantities (connected parts) scale differently with the volume than the parts of 
which they are made. For instance, for (3 < j3c and n sufficiently large. 



< >n -3(< M„ oc 2" , (8.1) 

while the individual components scale like 2^". The situation is worse if we consider the 
6-point functions, where the connected part has the same scaling as Eq. ( ^.Ij ) but the 
individual components scale like 2^". In other words, the beginning significant numbers of 
the individual terms do not matter for the subtracted quantities. Assuming 16 significant 
digits, when 2^" reaches 10^^, we still get the subtracted parts with 8 significant digits. 
When 2"" reaches 10^^, there are no significant digits left for the subtracted part. 

As a consequence, it is not always possible to stabilize the value of the connected part 
during as many iterations as we would like, given the study of section ^ This is an interesting 
situation. As long as we increase the number of iterations, we get a value of the unsubtracted 
quantity which becomes closer to its infinite volume limit. If we represent the significant 
digits of a double precision number as a sequence of 16 digits written in the conventional 
way, we can visualize this procedure as the successive obtention of the digits on the right side 
of the number. Unfortunately, at the same time, the part which gets subtracted increases 
in magnitude. Consequently, more an more digits on the left side of the number are wasted 
for the evaluation of the subtracted quantities. The situation gets worse if we consider the 
6 or higher point functions. 

The subtracted quantities are diverging near criticality. However, it is possible to define 
dimensionless renormalized coupling constants which have a finite limit. In the case of 
the four point function, the dimensionless renormalized coupling constant A4 is obtained by 
multiplying the zero-momentum connected Green's function 6*4 by the D + A power of the 
renormalized mass m^, namely 

A4 = -G'>g+^ . (8.2) 
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The mnemonic for D + 4 is 8 (amputation of the 4 legs at zero-momentum) + D — 4 (the 
canonical dimension of the (f)'^ bare coupling constant). We are using the notation 

G;.K„„_l<2^^1li^^ (8.3) 

In order to compare with field theory results, one should consider Landau- Ginzburg 
measures where the cut-off dependence has been restored explicitly. For instance, in D = 3, 
the definition of the functions entering in the Callan-Symanzik equations (the beta function 
etc.) given in Ref. [|r^ requires that we keep the dimensionful constant fixed while the cut- 
off goes to infinity. In other words, we need to change the dimensionless constant entering 
in Ro{k) while taking the infinite cut-off limit. This delicate procedure is beyond the scope 
of this paper, where we emphasize the basic numerical aspect of a single calculation. As 
explained in the introduction, we continue using a fixed Ising measure and a single adjustable 
parameter [(5). 

The quantity A4 has a finite (and supposedly non-zero when D < 4) limit when A — 00 
or equivalently when (3 — >• j3c. We can thus bypass the explicit introduction of the cut-off. 
Taking into account that there is no wave function renormalization, or in other words that 
the critical exponent rj is zero, we define A4 as the limit where n goes to 00 of 

_ < >„ -3(< Ml >„)^ 



Equivalently, with the convention of Eq. (|2.9| ), which does not involve inverse factorials, 
and for the value s = 2c~5 always used here, we obtain 

^'n 1 — '^^n 2 

In practice, we pick a given relative precision A and we require that n is large enough to 
stabilize the susceptibility and A4 with a relative precision 6. The reason for requiring both 
conditions is that A4 may temporarily stabilize when the flow passes near the fixed point 
(and so we are still far away from the infinite volume limit), but this is signaled by the fact 
that the susceptibility is still growing. In summary, we require 

.«"+i'i_^|<A (8.6) 



and 



O-nA 



4-"+i_i|<A. (8.7) 



A 



4,n 

When these two conditions are satisfied, we check that our result for A4 is compatible with 
the expected precision, or in other words, that we have enough significant digits left in 
I — 2a„,2 to calculate A4 with a relative precision A. We thus require the additional 
condition 
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where (5 is a typical round-off error (10~^^ in double precision). If the additional condition 
is not satisfied, we lower 6 and repeat the calculation. We have applied this algorithm in 
D = 3, 4, and 5 and for —Logio{Pc — /3) = 2, 3, . . . , 14. 

For D = 3, we were able to do all the calculations with A = 10^®. We found that A4 
reaches a limit A4 = 1.92786 when (3 —>■ (3c- In other words, hyperscaling holds very well. 
Fig. |T3| shows that to a good approximation 

A4 - a: ~ 1.68 X (8.9) 

In D = 4, we had to reduce to A = 10""^. We found that A4 tends to zero when (3 (3^- 
As shown in Fig. we have the approximate law 

1 



-1.96 - 0.746 X Ln{l3c - (3) 



i.lO) 



which is consistent with perturbative calculations JT^|. 

In D = 5, we had to reduce further to A = 10^^. We found that A4 tends to zero 
according to the approximate law 

A4 ~ 1.02 X (/3c - /5)°-^°^ (8.11) 



as shown in Fig. |T5|. If we replace {jSc — (3) hj A ^, we see that our result is consistent with 
A4 oc A"^ 



IX. CONCLUSIONS AND PERSPECTIVES 

We have shown that the use of truncations in the Fourier transform of the recursion 
relation of Dyson's hierarchical model leads to systematic errors which can be suppressed 
more than exponentially when the dimension of the truncated space increases. We have 



justified the use of the truncation for calculations |T0[ of the high-temperature expansion. We 
have shown that the finite volume effects can be reduced with an exponential precision. We 
have found the temperature dependence of the numerical errors and explained the empirical 
results with a simple renormalization group argument. 

The numerical errors appear as a practical aspect of the hierarchy problem. If it seems 
hard to believe that nature would fine-tune its fundamental parameters to produce scalar 
particles with masses very small compared to the Planck scale, there are also practical 
difficulties related to this fine-tuning. In the present context, it is quite simple to find 
f3c- However, this is not the end of the story: the physical quantities become numerically 
unstable when we reach j3c- This difficulty is not unsurmountable if we want to reach a 
cut-off of the order of the Planck scale which is only 17 orders of magnitude larger than the 
weak scale. We can use programming methods with enough significant digits. To take the 
analogy with differential equations, the problem of sensitive dependence on initial conditions 
can be dealt with provided that we do not evolve the system during a too-long amount of 
time. 

The methods presented here can be applied to field theoretical calculations. The simplest 
one is the calculation of the renormalized mass. This quantity is crucial because it enters 



14 



into the definition of tfie functions appearing in the Callan-Symanzik equations |T^. For 
the model considered here, a possible definition of the renormalized mass is 



mW, = , (9.1) 

where A is the largest eigenvalue of the linearized renormalization group transformation 
which needs to be calculated precisely. A is a UV cut-off taking the value 2"d Ar where Ar 
is a scale of reference below which we are considering an effective theory. Finally, yU is a 
parameter which allows us to change the value of the renormalized mass. As explained in 



section |VIII| , the method can also be applied to the calculation of renormalized quantities 
in Landau-Ginzburg models. These calculations will be used to check the validity of the 
perturbative evaluations of the functions entering the Callan-Symanzik equations. In the 



case D = A, the high-temperature expansion [|T^ indicates that the perturbative result |T5 
is very accurate. 

More generally, the calculational method presented here can be used to check any kind 
of approximate calculation which applies to the hierarchical model. 
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FIG. 2. -J- versus the dimension D. 
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FIG. 3. Relatives differences between the susceptibility calculated with the main method (x) 
and three alternative methods (x^*'') with i = 1 (crosses), 2 (diamonds) and 3 (circles) as in the 
text. The calculations were done in D = 3 and with (5 = f3c — lO"^-*^. 



19 




5 10 15 

-Log,o(/?^-/S) 



FIG. 4. Relative difference between the susceptibility calculated with the main method x 
with a rescaling x*^"^^ as explained in the text versus /3c — /3, in L> = 3 (circles), D = A (stars) and 
D = 5 (squares). 
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FIG. 5. n(/3, A) defined in Eq.( |5.2D , for D = 3 (circles), D = 4 (stars) and D = 5 (squares). 
Filled symbols correspond to a fixed value A = 10~^^, empty symbols correspond to a variable 
value A = IQ-^^ / {/3c - P) ■ 
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FIG. 6. Relative difference between tlic susceptibility calculated with Imax = ^ + 1 and Imax = I 
in D = 3 (circles), D = A (stars) and D = 5 (squares). /? = /3c — 10~^ in the three cases. 
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FIG. 7. Relative difference between tlic susceptibility calculated with Imax = ^ + 1 and Imax 
D = 3 (circles) and with /3 = /3c — 10~^, compared with a linear fit of these points. 
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FIG. 8. Relative difference between the susceptibility calculated with Imax = / + 1 and 
Imax = / in D = 3 with (3 = (3c- 10"^ (filled circles), P = Pc - 10"^ (asterisques), P = Pc - 10"'' 
(crosses), up to P = Pc — lO"-*^^ (empty circles). 
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FIG. 9. Intercept of the linear fits corresponding to the hnear part of Fig. 
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FIG. 10. Relative difference between tlie coefficients of tlie susceptibility b„i calculated with 
Imax = ^ + 1 and Imax = / in D = 3 with m=200 (squares), m=300 (circles) and m=400 (triangles). 
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m — dependence of s (D = 3) 
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FIG. 11. Ln{—s{m))/Ln{m) versus Ln{m) in D = 3. 
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m — dependence of s (D = 4) 




FIG. 12. Ln{—s{m)) / Ln{m) versus m in = 4. 
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FIG. 13. Logio{X4 — A4) versus —Logio{Pc — /?) in 3 dimensions. 
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FIG. 14. I/A4 versus —Logio{(3c — /?) in 4 dimensions. 
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FIG. 15. Logio{\4) versus —Logio{l3c — /?) in 5 dimensions. 
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